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ABSTRACT 

The bursting pulsar , GRO J1744-28, went again in outburst after ~18 years of quies¬ 
cence in mid-January 2014. We studied the broad-band, persistent, X-ray spectrum 
using X-ray data from a XMM-Newton observation, performed almost at the peak 
of the outburst, and from a close INTEGRAL observation, performed 3 days later, 
thus covering the 1.3-70.0 keV band. The spectrum shows a complex continuum shape 
that cannot be modelled with standard high-mass X-ray pulsar models, nor by two- 
components models. We observe broadband and peaked residuals from 4 to 15 keV, 
and we propose a self-consistent interpretation of these residuals, assuming they are 
produced by cyclotron absorption features and by a moderately smeared, highly ion¬ 
ized, reflection component. We identify the cyclotron fundamental at ~4.7 keV, with 
hints for two possible harmonics at ~ 10.4 keV and ~ 15.8 keV. The position of the 
cyclotron fundamental allows an estimate for the pulsar magnetic field of (5.27 ± 0.06) 
x 10 11 G, if the feature is produced at its surface. From the dynamical and relativ¬ 
istic smearing of the disk reflected component, we obtain a lower limit estimate for 
the truncated accretion disk inner radius, (>100 R s ), and for the inclination angle 
(18°-48°). We also detect the presence of a softer thermal component, that we associ¬ 
ate with the emission from an accretion disk truncated at a distance from the pulsar 
of 50-115 R 9 . From these estimates, we derive the magneto-spheric radius for disk 
accretion to be ~0.2 times the classical Alfven radius for radial accretion. 

Key words: line: identification - line: formation - stars: individual (GRO J1744-28) 
— X-rays: binaries — X-rays: general 


1 INTRODUCTION 

Mildly-recycled, slowly spinning, accreting X-ray pulsars in 
low-mass X-ray binaries (LMXBs) are uncommon objects, 
because of the relatively short time required to spin up the 
neutron star (NS) to frequencies f sp in > 100 Hz during the 
X-ray active phase. Before, or simultaneously, with this re¬ 
cycling phase, the magnetic field (B-field) of the NS decays 
from values B > 10 12 G, typical of young NS in high-mass 
X-ray binaries (HMXBs), down to B < 10 s G, when the 
strength of the field is no longer able to significantly drive 
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the accretion flow and coherent pulsations in the persistent 
emission are no longer detected. 

Our current understanding of h ow NS B-fields are actu¬ 
ally dissipated is not y et complete dPavne fe Melatosll2004l ; 
IZhang fe Koiimall2006l) . so that the study of LMXBs where 
measures of the B-field are possible are of extreme import¬ 
ance. Direct and indirect methods for the determination of 
the B-field are possible if pulsations are detected, as in the 
case of recycled LMXBs that host fast spinning NSs (the 
class of the accreting millisecond X-ray pulsars, AMXPs) 
and for a smaller group of LMXBs where the process of re¬ 
cycling is probably caught in its very early stage, and the 
NS spin frequency is much below 100 Hz (IPatruno et al.l 
120121) . Members of this class are e.g. the 1.7 Hz pulsar in 
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X1822-371, where the mass transfer is probably highly non- 
conservative and at su per-Eddington rates feurderi et al.l 
l201(ll : Iaria et ah| 1 201 3|) : the 11 Hz pulsar in IGR J1748- 


2446 (Pagitt^ et_alj 201111 : the 0.13 Hz pulsar in 4U 1626- 
67 dBeri et al. 2014h : the 0.8 Hz pulsar in the well-known 
system Her X-l ( Tananbamnjh i _al] jlj)72 1 and the 2.14 Hz 
pulsar in GRO J1744-28 llDeeenaar et abl2014h . 

GRO J1744-28 was discovered in hard X-rays (25- 
60 keV) by BATSE/C GRO on 2nd, December 1995 
dKouveliotou et al.lll996l l as a bursting source. Bursts were 
soon associated to unstable/spasmodic accretion episode s 
(type-II burs ts. Kouveliotou et al .11 1 9961 : Ike win et ahlll996lh 
iFinger et al.l (1 1996T) discovered coherent pulsations at 2.14 
Hz and its Doppler modulation, leading to the determ¬ 
ination of the orbital period (11.8337 days), eccentricity 
(e < 1.1 x 10 —3 ), projected semi-major axis (2.6324 lt-s), 
and mass function (1.3638 x 10 -4 Mg). The pulsed profile 
showed a nearly sinusoidal shape, which together with the 
value of the mass function, implies the inclination angle, or 
the companion mass, or both these conditions, to have low 
values. 


We observed to date three major X-ray outbursts from 
GRO J1744-28: the first associated with its discovery last¬ 
ing from December 1995 to May 1996, wher e at t h e peak 
outburst reached over-Eddington luminosities llJahoda et al.l 
Il999l : iGiles et ahlll996l : IWoods et al .111 999l . I 2 OO 0 I) : the second 
outburst occurred about one year later, started at the be¬ 
ginning of December 1996 and lasted about four months. A 
long period of quiescence followed, interrupted by the third 
and most recent outburst observed in January 2014. 


It was soon realized, just after its discovery, that the 
low spin period, the long orbital period and the observed 
flux and spin derivative are all factors that point to an in¬ 
termediate B-fiel d pulsar (10 11 G < B < 10 12 G) in a low in¬ 
clination system llDaumerie et aljfl996l 'l . ICuil (Il997l ~) presen¬ 
ted observational evidence for a propeller effect during the 
decay period of the first outburst, where pulsations were not 
detected below a flux of 2.3 x 10~ 9 erg cm~ 2 s -1 (2-60 keV), 
from which a rough limit on the Alfven radius could be de¬ 
rived, and thus leading to a B-fiel d estimate of ~ 2.4 x 10 11 
G, assuming a distance of 8 kpc. iRappaport fe Jossl (ll997T ) 
presented possible evolutionary paths for GRO J1744-28 
that could lead to the formation of the presently observed 
system, and starting from a set of reasonable initial con¬ 
ditions, gave rather tight limits on the possible companion 
mass (0.2-0.7 M@), neutron star mass (1.4-2.0 Mg), inclin¬ 
ation angle (7°-22°) and pulsar B-field (1.8-7 x 10 11 G). 


Search for the optical/NIR counterpart has been dif¬ 
ficult because of the crowded Galactic Centre field where 
GRO J1744-28 resides. Chandra and XMM-Newton de- 
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activity in 2002 dWiinands fe Wand 120021 : iDaigne et al.l 
2002), allo wing to refine t he est imates on the source co¬ 
ordinates. D egenaar et al.l (l2012l l reported enhanced flux 
Lx ~ 1.9 x 10 34 erg s _1 ) in September 2008. IGosling et al.l 
2007) determined two near-infrared possible candidate 
sources withi n the Chandra err or circle. Most recent de¬ 
terminations dMasetti et al. 2014l~) indi cate that the candid¬ 
ate source a by IGosling et alT(l2007f ) is the likely coun¬ 
terpart. Hence, the companion star is an evolved G/K III 
star. GRO J1744-28 has been associated to the Galactic 


Centre populatio n, _ setting its distance at ~ 7.5 kpc 

dAugusteiin et al1ll997h . 

One of the main peculia rity of GRO J174 4 -28 is the 
presen ce of Type-II bursts dLewin et aid 1 1 9961 1 . ICannizzol 
IU9971 1 identified the Lightman-Eardley instability occurring 
at the magnetospheric radius in the ac cretion disk as the 
possib le cause of the bursting behaviour. iBildsten fe Brown] 
d 199711 argued for the possible presence of type-I bursts 
together with the type-II as observed also for the Rapid 
Burster for accretion rates lower than 6 x 10~ 9 M@ yW 1 . 
The bursts have similar structures, with a fast rising (few 
seconds), a peak emission that lasts few seconds and a tail 
that appears as a dip with respect to the pre-burst average 
count rate. During the bursts the pulsed fraction increases 
but the phase with the pre-outburst pulsations is not main- 
tained, leading to the unusual phenomenon of pulsar glitches 
dStark et al.lll996l l. 

The X-ray continuum emission of GRO J1744-28 during 
the 1996 outburst has been st udied using the lo w energy- 
resolution RXTE /PCA data. I Giles et all dl996l l modelled 
this spectrum using an absorbed power-law with a high- 
energy cut-off (E cut = 20 keV, E/ 0 (d = 15 keV). They also 
found a pulse-dependent photon-index between 1.20 ±0.02 
(at the peak of the pulsed emission) and 1.35 ±0.01 (at 
the pulse minimum), and detected a bro ad (a ~ 1 keV) 
Gauss ian emission line in the Fe Ka region. iNishiuchi et al.l 
d 19991) presented a spectroscopic investigation with ASCA 
data for two different dates and accretion rates. They found 
the line to be intrinsically broad and they proposed a relativ- 
istically smeared reflection component as the cause of the 
broadening. The iron line profile was compatible with res¬ 
onant emission from mildly ionized iron for the first observa¬ 
tion, whereas highly ionized (6.7-6.97 keV) iron seemed pre¬ 
ferred for the second observation. However, the best-fitting 
value for the inclination angle, > 50°, seems unlikely be¬ 
cause of the constraints on the mass function and compan¬ 
ion’s spectral type. Similarly, the value for the inner disk 
radius, 10 R 9 < R™ < 30 Rjj) appears also unlikely because 
of the different indica t ions for an intermediate NS B-fiel d 
llDaumerie et al.l 1 19961 : ICuil 1 19971 : IRappaport fe Jossl Il997f ) , 
that would truncate the disk at larger distances, so that 
other broadening mechanisms, or a more complex profile, 
hav e been ad ditionally advocated. 

iDegenaar et al.l (120141 ') reported from a recent high- 
resolution Chandra/HETGS observation the detection of an 
asymmetric line profile from He-like iron, and interpreted it 
as a disk reflected feature, from which an inner disk radius of 
85 ± 11 R a an d an inc lination angle of 52° ± 4° have been es- 
timated. lYounes et al.l (l2015l l analysed the broadband spec¬ 
trum of GRO J1744-28 with NuSTAR and Chandra , fitting 
the continuum emission with a cut-off power-law of photon- 
index r ~0, cut-off energy ~7 keV, and a soft continuum 
thermal component at a temperature of 0.55 keV, inter¬ 
preted as emission from the accretion disk. They also repor¬ 
ted the presence of a significantly broadened (cr = 3.7 keV) 
feature at ~ 10 keV. Moderately broadened highly-ionized 
emission lines at 6.6 keV and 7.01 keV were also detected, 
together with emission from neutral iron at 6.4 keV. 


1 The gravitational radius R g = GM/c 2 is ~ 2 km for a canonical 
1.4 M 0 NS 
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Figure 1 . Swift -BAT jKrirnrn et al .1 [20~K?1 daily light curve in 
the 15-50 keV range from 22nd January to 1st June 2014. The 
black thick (dotted) arrow marks the time of the XMM-Newton 
{INTEGRAL) observation. 


The 2014 outburst of GRO J1744-28 

Enhanced X-ray emission from the direction of the transient 
outbursting source GRO J1 744-28 was detecte d from 18th, 
January 2014 with MAXI (jNegoro et al.ll2014l ). Pulsations 
at 2.14 Hz were detected in hard X-rays with F ERMI/GBM 
monitor on 24th, January 2014 (iFinger et al.l 120141 1 and in 
soft X-rays with the Swift/X RT on 6th, February 2014 
(iD’Ai et al.ll2014ll . confirming the association of the ongoing 
outburst with the bursting pulsar. The outburst reached its 
peak in mid-March, when the source was accreting close to 
the Eddington limit dNegoro et al.l [201411 and from then it 
started to decline linearly until 23th April when the light 
curve as seen in the Swift /BAT monitor (Fig. [T} clearly 
showed a drop, and the profile became exponential with a 
decay time of ~ 11 days. The outburst was r egularly mon - 
itored by Swift and INTEGRAL satellites dMasetti et ali 
l2014|j . A Target of Opportunity observation with XMM- 
Newton observed the source from 6th 14:03:12.27 to 7th 
13:29:50 March 2014, for a total exposure time of 83915 s. 
The observation was performed few days before the GRO 
J1744-28 flux reached the outburst peak as shown by the 
black thick arrow superimposed to the BAT light curve from 
MJD 56679 (22nd, January 2014) to MJD 56809 (1st, June 
2014) in Fig.|U For this work we also exploited data from an 
INTEGRAL observation (dotted arrow in Fig. [T]) performed 
three days later the XMM-Newton pointing. 
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instruments (RGS1 and RGS2) in Spectroscopy Mode. The 
EPIC/MOS instruments were switched off to allow the 
highest possible telemetry for the EPIC/pn. 

The RGS data were extracted using the rgsproc 
pipeline task. The total exposure for the RGS1 (RGS2) 
instrument is 83715 s (83750), and the first-order net 
light curves has an average count rate per second (cps) 
of 0.283 ±0.002 (0.268 ± 0.002). Second-order count rates, 
for RGS1 and RGS2, are ~30% and ~70% the first-order 
rates, respectively. The bursts are not clearly resolved in the 
RGS light curves, because the source is strongly absorbed 
at lower energies and the RGS1 and RGS2 frame times are 
4.59 s and 9.05 s, respectively, that are long compared to 
the rapid burst evolution. 

The EPIC/pn event file was processed us¬ 
ing the epproc pipeline processing task with the 
runepreject=yes, withxrlcorrection=yes, runepfast=no. 
and withrdpha=yes, as suggested by the most recen t 
calibration status repoiij;] (see also IPintore et al.l |2014| ). 
The average count rate, corrected for telemetry gaps 
(epiclccorr tool), during the entire observation over all 
the EPIC/pn CCD was 714 cps. The EPIC/pn light curve 
shows 43 bursts during the entire observation with a mean 
rate of 0.53 burst hr -1 . Fig. [5] shows the first 20 ks of the 
EPIC/pn observation, where nine bursts are present. The 
bursts profiles can be clearly resolved, they have similar 
structures with an average rising time and typical peak 
duration of a few seconds. During these peaks the observed 
rate becomes so high that telemetry saturation and strong 
pile-up prevent a reliable analysis of the data. 

The burst decay is well described by an exponential tail 
with decay-times in the range 2-10 s. The count rate after 
each burst drops about 10% with respect to the pre-burst 
level and returns to it slowly after few hundred seconds. Such 
behaviour and the burst characteristics are very close to the 
typical ph enomenology show n by the source in the two past 
outbursts dGiles et al. Il996l 'l. 

In Fig. [3] we show a snapshot of the EPIC/pn light 
curve encompassing two typical bursts (upper panel), and 
the hardness ratio (HR) evolution, where the soft color is 
defined in the energy range 0.5-4 keV and the hard color 
in the 4-10 keV range (lower panel). The light curve is cor¬ 
rected for telemetry saturation (epicclcorr task), and it is 
filtered removing the three hottest RAWX columns to limit 
pile-up effects at the burst peak. 

We focused on the persistent emission of the pulsar by 
time-filtering the EPIC/pn event list excluding appropriate 
intervals around the bursts. The excluded windows are all of 
400 s length, starting a few seconds before the rising trail of 
the bursts. The final exposure, corrected also for telemetry 
drops, of the persistent emission resulted of 64130 s. 


2 OBSERVATION AND DATA REDUCTION 

We used the Science Analysis System (SAS) version 13.5.0 
and the calibration files (CCF) updated at March 2014 for 
the data extraction and analysis of the XMM-Newton data; 
we used the software package HeaSOFT version 6.15.1 for 
the timing analysis and for general purposes data reduction, 
and Xspec v. 12.8.1 for spectral fitting. 

The EPIC/pn operated in Timing Mode with the op¬ 
tical thick filter, while the Reflection Grating Spectrometer 


2.1 Timing analysis 

We barycentred the EPIC/pn, burst-filtered, event file with 
respect to the Barycentric Dynamical Time (TDB) using the 
barycen tool with R.A. and DEC coordinates 266.137792 
and -28.740803, respectively ll Gosling et al]l2007l) . We used 

2 http://xmm2.esac.esa.int/docs/documents/CAL-TN-0018. 
pdf 
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Figure 2. EPIC/pn light curve for the first 20 ks of the observa¬ 
tion. Bin time is 5 s. 
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Figure 3. Two consecutive bursts and relative hardness ratio 
(HR, 4.0-10.0/0.5-4.0 keV). Selected extraction region excludes 
the RAWX= [36:38] columns to avoid photon pile-up. Bin time 
is 5 s. 
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all the events that passed the filtering condition FLAG==0. A 
power spectrum of the persistent emission clearly reveals 
the pulsation at 2.14 Hz and the first harmonic at 4.28 
Hz. To correct the pulse arrival times for the orbital mo¬ 
tion of the pulsar, we used the following orbital parameters 
(GBM pulsar projec10): P or b = 11.836397 d, a sin(i) = 2.637 
lt-s, and Tt /2 = 56695.6988 MJD. We then used a folding 
search (efsearch in XRONOS) to find the averaged pulse 
period during the persistent emission. The \ 2 highest peak 
is found in correspondence with a period of 0.4670450(1) s, 
or 2.14112130(46) Hz, where the uncertainty corresponds to 
the step period search time (10~ 7 s). 

To check possible shifts in the spin frequency during the 
observation, we folded at this period the persistent event file 


3 http: //gaimnaray.msfc.nasa.gov/gbm/science/pulsars/ 
lightcurves/groj1744.html 



Phase 

Figure 4. Top panel: energy-averaged pulse-folded profile (32 
channels) and best-fitting model using two sinusoids. Bottom 
panel: pulse-folded profile for the EPIC/pn data in energy bands: 
0.5-3.0 keV (red data), 3.0-6.0 keV (green data), and 6.0-12-0 
keV (blue data). 


using a folding time window of 10 4 pulse periods and the 
start time of the observation as epoch of reference. We then 
plotted the phases with respect to this folding period as a 
function of the time. The phases showed random scattering, 
and no other correction was found necessary. Folding the en¬ 
tire persistent emission at this period resulted in the folded 
profile of Fig. [I] The pulse shape can be well fitted as the 
sum of two sinusoids. Using the standard notation for the 
pulse amplitude A, 

. I max I min 

A ~~ 7 _L T ■ 

-t max i J-rmn 

we found the fundamental has an amplitude of 15.1%, while 
the first harmonic amplitude is 0.91%, corresponding to 
~ 6% the fundamental. The folded profile is strongly energy 
dependent in the EPIC/pn band, as clearly shown in the 
bottom panel of Fig. [I] where the folded profile is obtained 
for three energy band: 0.5-3.0 keV, 3.0-6.0 keV and 6.0-12.0 
keV. 
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Figure 5. Amplitude fraction for the fundamental (red data, left 
y-axis) and first harmonic (green data, right y-axis) as a function 
of the energy. 


Figure 6. Upper panel: fundamental amplitude as a function of 
energy in the 2.0-12.0 keV range, super-imposed with the best¬ 
fitting trend derived by data in the same range, excluding the 
interval 5.5—7.5 keV. Lower panel: residuals in units of sigma for 
the 5.5-7.5 keV range for the best-fitting model consisting of a 
third degree polynomial derived by fitting the data in the 2.0—5.5 
keV and 7.5—12 keV range. 


2.2 The amplitude fraction 

We studied the amplitude fraction as a function of energy 
for the EPIC/pn band. We exploited all the selected events 
in the RAWX columns, although we were aware that pile-up 
may distort the energy-dependent amplitude measures; how¬ 
ever we estimated that the relative error induced by some 
fraction of pile-up events is significantly lower than the stat¬ 
istical uncertainties introduced when removing the hottest 
columns. We selected the pulse profiles in the pulse-invariant 
(PI) channel as follows: from 0.5 to 2 keV by a step of 0.5 
keV, from 2 keV to 8 keV by a step of 0.25 keV, from 8 to 
10.5 keV by a step of 0.5 keV, and from 10.5 keV to 12 keV 
into a single profile. 

We folded the persistent event file using our best period 
estimate (0.4670450 s) and fitted the resulting pulse shape 
using two sinusoids. In Fig. [5] we show the dependence of 
the amplitude of the fundamental and first harmonic as a 
function of the energy. 

It is evident that above 2 keV the pulse amplitude 
strongly correlates with energy, with the exception of the 
energy range between 6 and 7 keV where a remarkable drop 
in amplitude is present. To better characterize the shape of 
this drop we used a phenomenological model to describe the 
amplitude dependence in the 2.0-5.5 and 7.5-12 keV range. 
We described this trend using a polynomial fit, attaining a 
X 2 ~ 1 with a third order polynomial. We show in Fig. [6l 
the residuals of this fit when the data of the 5.5-7.5 keV 
region are then noticed. The shape of the residuals is con¬ 
sistent with a Gaussian profile, that it is centred at ~ 6.5 
keV and has a sigma of ~ 0.8 keV. This drop indicates the 
presence of a spectral component that is not pulsed, and we 
shall show that its shape is compatible with the shape of 
a disk reflection component as derived from the analysis of 
the broadband spectrum (see Sect. 13.51 . 


3 SPECTRAL ANALYSIS 
3.1 RGS data analysis 

RGS count rates are always below the threshold for photon 
pile-up. Because of the higher spectral resolution of the RGS 
detector with respect to the EPIC/pn and the low-energy 
coverage in a region where the effect of ISM absorption is 
the highest and the EPIC/pn spectrum is affected by larger 
systematics, we first performed a preliminary spectral ana¬ 
lysis in this restricted energy band. We extracted first and 
second order RGS spectra of the persistent emission using 
the same good time intervals obtained from the EPIC/pn 
light curve. We looked at the background contribution to 
choose the noticed energy ranges for spectral fitting, keeping 
channels above ~ 3 a the background value. Using a simple 
model consisting of an absorbed power-law and a constant 
of normalization free to vary for the different datasets, we 
found this inter-calibration factor among the four spectra 
variable at the 7% level. 

We verified the consistency of the RGS1/RGS2 spec¬ 
tra and generated a combined spectrum using the SAS 
rgscombine tool for the first (RGSol) and second (RGSo2) 
order spectra. After having verified that the persistent spec¬ 
trum is, within the statistical uncertainties, consistent with 
the total time-averaged spectrum, to increase the statist¬ 
ical weight of the RGS data, we used this time averaged 
spectrum for the rest of our analysis. We found no evid¬ 
ent emission/absorption narrow feature, so that we chose 
to group the spectra to a minimum of 100 counts/channel. 
Both spectra (first and second orders summed) are used in 
the 1.3-2.0 keV range, because below 1.3 keV the spectrum 
becomes background-dominated. 
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3.2 EPIC/pn data analysis: pile-up and 
background issues 

Because of the high count rate recorded by the EPIC/pn 
instrument, we checked the pile-up level of the persistent 
time-averaged emission. In the EPIC/pn Timing mode spa¬ 
tial information is collapsed into one dimension (RAWX co¬ 
ordinate), while spectral extraction is usually performed in 
a 2-dimensional array, using the RAWY coordinate, that, 
however, is a temporal coordinate. 

The point-spread function in the RAWX-RAWY im¬ 
age peaks for the RAWX=37 column and we created four 
filtered spectra selecting the RAWX=[28:46] region, and suc¬ 
cessively removing one, three, and five central columns, ac¬ 
companied always by the filtering expressions PATTERNS 
and (FLAG==0) . We obtained response matrices using the 
rmfgen tool and the ancillary files using the art gen tool, 
following standard pipeline^. 

For each filtered event list, we applied the epatplot tool 
and we noted that the exclusion of three central columns 
provided a spectrum with an observed-to-model singles and 
doubles pattern fraction close to unitjjf]. However, we noted 
that the choice of excising the PSF caused severe effects 
on the continuum spectral determination. In fact, when 
the spectra were comparatively examined, it resulted evid¬ 
ent that the removal of central columns softened the spec¬ 
trum beyond what could be expected by the pile-up cor¬ 
rection. We probed this by means of spectral analysis. In 
Fig[T] we show the spectra extracted from the EPIC/pn 
using all the RAWX[28:46] columns ( full spectrum ), and 
removing the three hottest columns, together with a ~ 1 
ks Swift/X RT spectrum (WT mode) taken on March 7th 
(OBS.ID 00030898022) in a simultaneous time window with 
XMM-Newton. We note here that the averaged count rate 
for the Swift/X. RT spectrum was ~70 cps, which is con¬ 
siderably below the threshold limit (~ 100 cps) for the WT 
mode photon pile-up, and the spectrum can be considered 
a bona fide benchmark. 

The excised spectrum clearly displays a strong flux drop 
above 7 keV that is not observed in the Swift/X RT spec¬ 
trum. The color ratio of the 2.0-6.0/6.0-10.0 keV absorbed 
flux is 0.61 for the spectrum without excision, and 0.69 for 
the spectrum with three columns removed, while for the 
Swift/XRT spectrum we measured a value of 0.63, that is 
much more consistent with the spectrum with no removed 
column. We argue that, although some level of pile-up could 
be present, the choice of excising the PSF artificially intro¬ 
duces spectral distortions, possibly related to the reconstruc¬ 
tion of the energy-dependent ancillary response, that be¬ 
comes more significant than the eventual pile-up correction. 
We choose therefore to use the whole PSF extraction region 
for spectral analysis, but we will also show the best-fitting 
results for the three-columns excised spectrum, to have a 
rough estimate of the sensitivity due to these instrumental 
issues on the determination of the physical parameters. Be¬ 
cause the EPIC/pn full spectrum contains ~5xl0' counts, 
statistical errors become possibly dominated by systematic 

4 http://xmm.esac.esa.int/external/xmm_user_support/ 
documentation/sas_usg/USG/epicpileuptiming.html 

5 http://xmm.esac.esa.int/sas/current/documentation/ 
threads/epatplot.shtml 
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Figure 7. Data, unfolded spectra and data/model ratio for the 
EPIC/pn spectrum extracted from the whole PSF (no columns 
removed, black data), for the EPIC/pn spectrum extracted from 
the wings of the PSF (three columns removed, red data), and 
for the Swift/XRT spectrum from a simultaneous observation 
(green data). 

errors of the response matrix and, to account for this, and 
obtain final y 2 close to 1, we adde d a 0.5 % systematic e rror 
to the EPIC/pn channels (see also iHiemstra et al.ll20lJ) . 

It is suggested to extract the background spectrum from 
the first columns (e.g. selecting RAWX columns between 2 
and 6) of the RAWX-RAWY image, however the flux of the 
source is severe, causing considerable spreading of the point 
spread function leading to contamination of source photons 
even in these first columns. Light curves extracted from this 
region and from central columns do in fact show a very sig¬ 
nificant positive correlation (p-value < 2.2e-16). To obtain a 
measure of the expected level of background, we used a 2001 
XMM-Newton observation covering the field of view of GRO 
J1744-28 (Obs.ID 0112971201) of ~ 3ks performed with the 
EPIC/pn in Imaging Mode (Full-Frame). In this observation 
the source is quiescent and no other point source is detected 
in a STegion around the source coordinates. We extracted a 
spectrum around the coordinates of the source using a cir¬ 
cular region of 2.4'radius, thus covering the same geometric 
area of our selection. 

We fitted this spectrum using a phenomenological 
model of a broken power-law, obtaining a satisfactory de¬ 
scription of the data with a reduced \ 2 — 1-1 for 112 d.o.f. 
We used this best-fitting model to generate a faked back¬ 
ground spectrum using the same response matrix, ancil¬ 
lary response, and time integration of the source spectrum. 
We thus verified that the expected background contribu¬ 
tion compared with the observed source flux is very small 
in all the EPIC/pn range, contributing ~ 1% below 2 keV 
and ~ 0.2% above. The background spectrum extracted from 
the first columns contributes ~ 2% above 2 keV, so although 
some variability in the region around the source can be ex¬ 
pected, we argue that the background subtraction may af¬ 
fect the source flux determination with this level of uncer¬ 
tainty. We finally verified that setting no background for the 
EPIC/pn spectrum had no significant impact for the determ¬ 
ination of the spectral parameters reported in the analysis. 
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Finally, we re-binned the spectrum in order to have at 
least 25 counts/channel and set the oversampling at 3 chan¬ 
nels/energy resolution with the specgroup tool. We quoted 
spectral errors at Ax 2 = 2.7, if not stated otherwise. 

3.3 Spectral Analysis of the Persistent 
Phase-Averaged Emission 

We first studied the persistent, pulse-averaged, spectrum. 
For spectral analysis, we used the 1.3-2.0 keV range of the 
RGS spectra and the 2.0-9.5 keV range for the EPIC/pn 
spectrum. We first noted that the EPIC/pn data below 
2 keV do not well match the RGS spectrum, and this 
mismatch can be modelled by an non-physical soft excess 
present in the EPIC/pn data that is not required in the 
RGS. Such mismatch has been reported in many other Tim¬ 
ing Mode observ ations of highly a bsorbed (Nh > 10 22 cm -2 ) 
sources (see e.g. [D’Ai et al.l l2Q10h . and it has been investig¬ 
ated in one of the most recent calibration technical reports of 
the EPIC/pn team (CAL-TN-0083). Although we carefully 
applied the suggested corrections to the EPIC/pn processing 
pipeline, still the residuals do clearly point to the presence 
of this instrumental artefact, probably enhanced by the high 
statistics of the EPIC/pn spectrum. Because the spectrum 
below 2 keV is well covered by the RGS data (Fig. O, that 
do not suffer from any known instrumental issue, we chose 
to discard EPIC/pn data below 2 keV. Because the source is 
highly absorbed, we found that the choice of the abundance 
table is quite relevant for determining the absolute value of 
the equivalent hydrogen column and the slope of the con¬ 
tinuum emission as well, given the high correlation between 
these parameters. We used the tbabs model and tried all the 
different abundance tables av ailab le for Xsp ctfl setting the 
cross-section table to that of lYerner et al.l dl996i l. We de¬ 
rived first some estimates using only the RGS datasets and 
a simple power-law continuum, obtaining values for equi- 
val ent column between (6.24± 0.07) x 10 22 cm -2 , adopting 
the lAnders fe Grevessel (Il989l') tab l e, to (9.68 ± 0.40) x 10 22 
cm~ 2 . adopting the I Wilms et al.l J200d l abundance table. 
We then obtained the corresponding best-fitting values us¬ 
ing only the EPIC/pn band and finally chose to adopt the 
abundance table that provided the best match between the 
RGS and EPIC/pn energy bands for the rest o f our a nalysis 
(i.e. abundance aneb, from I Anders fe Ebiharal dl982l V). 

Because the spectrum of the source extends well above 
10 keV and the cut-off energy of the spectrum is close to the 
edge of the EPIC/pn band, it is extremely useful, to bet¬ 
ter constrain the overall X-ray continuum shape, to comple¬ 
ment the XMM-Newton data with observations from other 
satellites covering the hard X-ray emission of the source. 
To this aim we used publicly available INTEGRAL/lSGRl 
(20.0-70.0 keV range) and INTEGRAL/ JEMX (we used 
both instruments JEMX1 and JEMX2 in the 7.0-18 keV 
range) of the INTEGRAL observation from March 10th 
(Revolution #1392), three days after the XMM — Newton 
observation (see Fig. [l] where the Swift/ BAT rate is 
0.181 ±0.006 and 0.182 ±0.006 cts cm -2 s^ 1 for the dates 
of the XMM — Newton and INTEGRAL observations, 

6 https://heasarc.gsfc.nasa.gov/xanadu/xspec/manual/ 
XSabund.html 


respectively). We used a calibration constant to take into 
account flux variations occurred between the INTEGRAL 
and XMM-Newton observations, assuming that the spectral 
shape has not significantly varied. 

We started selecting the possible continua based on 
the physical characteristics of GRO J1744-28, that is quite 
likely an intermediate B-field source, accreting via Roche- 
lobe overflow, with presence of strong pulsed emission (15% 
in the average 1.0-10 keV range). The pulsed fraction emis¬ 
sion is intermediate between typical values found in HMXB 
pulsars and the ones found in AMXPs, and this supports 
the hypothesis that the NS magnetic field is not so high to 
efficiently capture matter from large distance in the disk. 
We expect therefore that part of accreting matter is chan¬ 
nelled and up-scattered in the accretion column, and that a 
(significant?) part of the X-ray emission may be produced 
either at the magnetospheric boundary or directly linked 
with emission from the accretion disk (if the disk is trun¬ 
cated relatively not far away) or from reflection by the ac¬ 
cretion disc of the primary radiation. 

However, if the magnetospheric radius R mo9 is 100 
R g , we expect a spectrum closer to that shown from HMXBs. 
Consequently, following the phenomenological descriptions 
mostly adopted to model the spectra of accreting HMXBs, 
we tried to fit the data with a set of commonly used spec- 
tral models for accr eting X-ray pulsars (see model details in 
ICoburn et al.l[20 02l: a cut-off power-law model (cutoffpl), 
a power-law with a high-energy cut-off model (highecut), a 
power-law with a Fermi-Dirac cut-off (f dco), and a negative¬ 
positive cut-off power-law model (npex). 

Alternatively, assuming the magnetic field is relatively 
low and unable to stop the accretion flow at large distance 
(Rmag < 100 Rg) from the NS, we fitted the data with two- 
component models as commonly done for low-inclination 
angle and high -M accreting LMXBs. 

We used for the hard X-ray emission from the accretion 
column a therm al Comptonization model, (nthcomp model 
IZvcki et alj 1 1999l l, and for the soft X-ray emission either 
a black-body component or a multi-coloured disk emission 
(ezdiskbb component in Xspec, model disk+nthcomp). 

We generally found very high reduced x 2 (>> 10), irre¬ 
spective of the continuum choice, with a common trend in 
the residuals, driven by the high signal-to-noise-ratio of the 
EPIC/pn channels. The residuals pattern has a clear broad 
peak centred in the 6.4-7.0 keV Ka iron range, a residuals 
dip from ~ 4 to ~ 6 keV, and some local peaks at ~ 2.2 keV, 
2.7 keV, 3.3 keV and 4.0 keV. The residuals of the iron peak 
strongly pivot the fit and cause a general softer continuum 
and a poor determination of the energy cut-off. We show in 
FigEl th e EPIC/pn residuals for three representative mod¬ 
els. 

Because the residuals cover a large energy range, it is 
not straightforward to determine their origin from a single 
corresponding local process. To simplify the presentation of 
our results, we show and discuss them in the frame of one 
representative continuum model; we chose the disk+nthcomp 
model because it provided the marginal final best x 2 an d it 
gives a reasonable representation of the main physical char¬ 
acteristics of the spectrum. We repeated the analysis using 
the set of continua listed above, and we found only small 
variations in the spectral parameters that did not affect the 
overall picture here presented. 
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Figure 8. EPIC/pn 2—9.5 keV residuals for three representative 
best-fitting continua. 


We first assumed that the major processes are in emis¬ 
sion and tried to model the residuals with a set of local emis¬ 
sion Gaussian lines, starting by the highest peak found in 
the iron range. Because the residuals indicated a broadened 
feature, and a single Gaussian proved inadequate to fit it, 
we assumed a possible blending of different transitions and 
we checked this hypothesis using two Gaussian lines with 
fixed energies at 6.7 keV (He-like iron) and at 6.97 keV (H- 
like iron), tying their widths, assuming they are produced 
in the same photo-ionised plasma. Both lines are strongly 
required, with a common width of ~ 0.35 keV (FigO panel 

A) , although the residuals landscape was still far from being 
satisfactory. 

We further examined the secondary pattern of local re¬ 
siduals outside the iron range. Because the peaked emission 
in the iron range had suggested presence of fluorescence lines 
of highly ionized elements, we followed this scenario and 
looked at the expected energies of H-like and He-like trans¬ 
itions of the most abundant elements (Si, S, Ar, Ca, and Fe). 
We added a set of Gaussian lines fixing the energies at the 
expected laboratory frame. To limit the degrees of freedom 
of the fit we set a common line width for the Si, S, Ar and 
Ca transitions and another one for the Fe transitions. Be¬ 
cause of the intensity of the iron lines, we also took into ac¬ 
count the possible presence of the K/3 transitions. We found 
significant detections (> 2 <r) for the following transitions: 
Sxvi (2.62 keV), Arxvm (3.32 keV), Caxix (3.90 keV), K a 
Fexxv (6.7 keV), Ka Fexxvi (6.97 keV), K/3 Fexxv (7.88 
keV) K/3 Fexxvi (8.25 keV). We also found a significant 
narrow line detection at ~ 2.2 keV This line is found in cor¬ 
respondence with a strong derivative of the effective area at 
an instrumental Au M-edge, and we retain it of instrumental 
origin, so that we will keep them for the evaluation of the 
fits but we will not longer discuss it. 

This best-fitting model and the residuals appeared, 
however, still unsatisfactory, because of the persistent S- 
shaped curvature between 4 and 6 keV (see Fig[9] panel 

B) . A possible physical process able to impress on the spec¬ 
trum such features is the cyclotron resonant scattering pro¬ 
cess (CRSFs). CRSFs, when besides the fundamental, higher 


harmonics are detected, are present in the spectra of young 
X-ray pulsars hosted in HMXBs, and are usually detected 
between 20 and 60 keV, whereas their detection in the spec¬ 
tra of X-ray pulsars in LMXBs is strongly hampered because 
of lower B-fields strengths and/or possible different accretion 
geometries. 

The position of the line, E cyc \ is related to the strength 
of the field by the simple relation: 

-Ecyd = 11.6 x B 12 /(l + z) keV (2) 

where the B-field is expressed in units of 10 12 G and the 
gravitational red-shift z is a function of the distance where 
the line is effectively produced (z ~0.3 for a line originating 
very close to the NS surface). For the case of GRO J1744-28, 
where a possible intermediate field between 1.8 and 7 x 10 11 
G has been claimed, the presence of a line between 1.8 and 
6.2 keV would be not surprising. 

To model this feature we used an absorption Gaussian 
line profile (gabs component in Xspec) that has three free 
parameters: the line energy position (Ei), the width (Wi) 
and the depth, Di, related to the optical depth, r, by the 
relation r = QAxD\/W\. We found a strong improvement in 
the x 2 f° r the addition of this component (the x 2 /degrees 
of freedom [dof] passing from 632/524 to 571/521), with an 
F-test probability of chance improvement of 2 x 10~ n . Fur¬ 
ther, we also found possible hints for the presence of higher 
harmonics at ~ 10 keV and at ~ 15 keV, however the depth 
and the position of these features showed strong correla¬ 
tion with other spectral parameters, resulting in large relat¬ 
ive errors. To limit these correlations and the implied large 
uncertainties, we imposed the width of the first (second) 
harmonic to have two (three) times the value of the corres¬ 
ponding fundamental width, and with this constraint, we 
found a reduced x 2 of 1.052 for 517 dof, resulting, how¬ 
ever, in a marginal detection at ~ 3.8 <r level. Finally, we 
noted a very narrow emission line, with an equivalent width 
of 9 eV, close to the energy of the neutral iron edge. The 
addition of this feature was found very significant (final re¬ 
duced x 2 1.021 for 515 dof, panel D of FigO. After this last 
component addition, no other evident feature was appar¬ 
ent in the residual pattern. The values of the instrumental 
constants for the different data-sets resulted generally reas¬ 
onable. There is a ~ 10% mismatch between the EPIC/pn 
data and the RGS first-order/second-order data, and it is 
most likely due to uncertainties in the reconstruction of the 
EPIC/pn effective area and to the EPIC/pn background is¬ 
sues. A cross-calibration factor, around 30%, is obtained 
between the EPIC/pn and the INTEGRAL/ JEMX2 and 
I NT EGRAL /ISGRI data, but these data are not simultan¬ 
eous with the EPIC/pn data as the INTEGRAL data were 
taken few days later the XMM — Newton observation. 

The continuum emission is characterized by three tem¬ 
peratures. We associated the softest (kTdisfc =0.54 keV) to 
the maximum temperature reached in a multi-coloured ac¬ 
cretion disk, and we derived from the normalization of the 
component (N ez ), an estimate for the accretion disk inner 
radius (R ez ) through the relation: 

Nez = cos(i) (-fTjz) ( 3 ) 

where D is the source distance and / is the color to 
effective ratio. We assumed for this estimate a D = 7.5 kpc, 
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i = 20° and / = 2, obtaining a best-estimate for R ez of 27.5 

Rg. 

The second temperature (kTw, = 1.8 keV) is the soft 
seed photon temperature and it is related to the radiation 
field which provided the source of soft photons, and it could 
correspond to the plasma temperature in the post-shock re- 
gion, or the mould t emperature at the surface of the NS 
iBecker &; Wolflll2007l ). The highest temperature (kT e = 6.7 
keV) is the thermal electron temperature, responsible for 
the Compton up-scattering of the soft photons (neglecting 
the contribution from bulk-motion Comptonization) and it 
is associated with the spectral cut-off energy. If we con¬ 
sider the RGSol flux as a more reliable proxy for estim¬ 
ating the source flux (because it is likely much less affected 
than the EPIC/pn spectrum by background and systematic 
uncertainties), we derived a 1.0-10.0 keV unabsorbed flux 
of 1.32 x 10“ 8 erg cm~ 2 s -1 , and an extrapolated 0.1-100 
keV flux of 3.0 x ICC 8 erg cm -2 s -1 . For a distance of 7.5 
kpc, the isotropic source luminosity is ^2.1 x 10 38 erg s -1 , 
slightly above the Eddington limit for a canonical 1.4 Mg 
NS (L Edd ~ 1.8 x 10 38 erg s -1 ). Within the energy band 
covered by the data (1.3-70 keV), the relative flux contribu¬ 
tion of the disk component is 6.2% of the total X-ray output 
from source. 

The final best-fitting model and residuals are shown in 
panel D of Fig. [9j the best-fitting parameters and errors 
of the continuum emission are shown in Tabic [T| and the 
best-fitting parameters and errors for the set of emission 
lines in Table [2] We determined the uncertainties on the line 
energies and widths using the EPIC/pn data set without any 
systematic error, because these parameters are not expected 
to be driven by this type of uncertainty. 

3.4 Modelling the highly-ionized features with a 
reflection component 

We have shown in the previous section as a possible inter¬ 
pretation for the large deviations from the continuum shape 
can be due to the presence of large, moderately deep, ab¬ 
sorption features that we interpreted as signs of CRSFs. We 
have also shown a complex pattern of local emission features, 
moderately broadened, and compatible with a set of highly 
ionized transitions from a photo-ionized environment. We 
investigate the characteristics of the line emission spectrum 
applying a self-consistent and physically motivated scenario 
for this set of local features. Excluding a possible origin in 
a thermal environment, that would require unrealistic MeV 
temperature, we make the hypothesis that this set of lines 
is due to reflection in an accretion disk that is still relat¬ 
ively close to the NS to produce appreciable dynamical and 
relativistic broadening. 

To this aim, we slightly changed the continuum shape in 
order to self-consistently incorporate the reflection compon- 
ent, adopting the reflection model re lxill (see details in 
iDauser et alll20fol ; lGarcfa et al. i[H3). This model assumes 
as incident continuum a cut-off power-law and it correctly 
takes into-account the angle-dependent reflection pattern. 
The reflection component is calculated only down to values 
of photon-index, T, of 1.4 and cut-off energy of 20 keV; such 
values are sensibly different from the continuum parameters 
that we found, however, because the reflection component is 
mainly a contributor in the iron range (fluorescence lines), 


Table 1 . Spectral fitting results for the model disk+nthcomp. 
Continuum and cyclotron line parameters. Errors given at Ay 2 = 
2.706 (90% confidence for the single interesting parameter). 


Parameter (Units) 

Value 

N h (10 22 cm" 2 ) 

6.3T0.1 

r 

1.96T0.09 

kT e (keV) 

6.75T0.26 

kT bb (keV) 

1.80T0.08 

Norm. a 

0 11+004 
U-1 - 0.01 

kT disfc (keV) 

0 54 "*”°' 11 
u - 04t — 0.07 

Inner disk radius R ez (R 5 ) 

27±7 8 

Ei (keV) 

c q 1 +0.09 
°- ol -0.14 

Wi (keV) 

0 39 + 0,18 
u.o»_o 09 

Di (x lO" 2 ) 

2 9 +0 ' 4 
z,y - 0.1 

E 2 (keV) 

11 4 +1 -° 

W 2 (keV) 

0.78 6 

D 2 (x 10" 2 ) 

22±10 

E 3 (keV) 

15 +20 
lo - 1.0 

W 3 (keV) 

1.17 c 

D 3 (x lO" 2 ) 

OZ -15 

Constant RGSol/PN 

1.13T0.02 

Constant RGSo2/PN 

1.08T0.02 

Constant JEMX1/PN 

1.18T0.02 

Constant JEMX2/PN 

1.32T0.02 

Constant ISGRI/PN 

1.33T0.10 

X 2 red ( dof ) 

1.00 (515) 


° nthcomp normalization, unity at 1 keV for a norm of 
b Parameter tied to 2 x Wi- 
c Parameter tied to 3 X Wi. 

we retain the model still a reasonable and satisfactory ap¬ 
proximation of the true reflection pattern. We left as free 
parameter the inner disk radius (R, r Jj) , the inclination angle 
and the emissivity index (e). The outer disk radius was held 
fixed at 1000 R g , and the spin parameter a=0. The other 
free parameters of this component are the ionization para¬ 
meter (log(£)), the reflection fraction ( Refl. fr.) and the 
iron abundance relative to the solar value (4i? e ). 

We tested our constraints on the widths of the CRSFs, 
but we found the widths of the first and second harmonic 
to be significantly different from 2 x Wi and 3 x Wi, and, 
to obtain a satisfactory fit, we left the CRSFs widths free 
to vary. We note that the two harmonics become statistic¬ 
ally more necessary with respect to the continuum model 
of Table Q] (both features detected at ~6.8cr, but the single 
harmonics are detected at ~3<r level), implying a certain 
model-dependent detection. We tested the presence of the 
softer disk component, left free to vary energy and normal¬ 
ization of the 2.22 keV instrumental feature. The best-fitting 
model has a final reduced chi-squared of 1.32 (526 dof), it 
describes quite satisfactorily the residual pattern in the iron 
range although we noted (left panel of Fie llOl) that some 
residuals still persisted at the energies corresponding to the 

7 To be noted that we drew another independent estimate for 
the same quantity from the best-fitting model derived from the 
normalization of the disk component, that we labelled R e z, Eq. [3] 
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Figure 9. Data, best-fit unfolded models (E X f(E) representation) and residuals in units of a for a model composed of an accretion 
disk component and a thermal Comptonization model. RGSol data in red, RGSo2 data in green, EPIC/pn (full spectrum ) data in black, 
JEMX1 data in blue, JEMX2 data in light blue, ISGRI data in magenta. RGS data graphically re-binned for clarity. Panel A: addition 
of two moderately broadened emission Gaussian lines in the iron range; panel B: addition of a set of highly ionized emission lines; panel 
C: addition of a gabs at ~ 5 keV. panel D: addition of two other gabs harmonics at ~ 10 and ~ 15 keV and of a zero-width line at 7.11 
keV. 


K a Si xvi, Ar xvm, Caxix and Caxx emission lines. To 
model the larger intensities suggested by these residuals, 
we used blurred relativistic profiles (relconv model) for the 
non-iron highly ionized emission lines, setting the blurring 
parameters (inner and outer disk radius, inclination angle 
and emissivity profile index of the relconv component) to 
be tied to the parameters of the relxill model. We found 
that the broadened line at ~3.9 keV, identified with the 
Ka Caxix transition in Table 0 is much better modelled 
with a combination of He- and H-like calcium emission lines. 
We found that the energies of the Ka Arxvm, Caxix and 
Caxx are consistent with the expected rest-frame values, 
but the position of the S xvi line is found at higher energies 
(2.73 ± 0.02 keV), however setting the line energy at 2.62 
keV still gives a lower limit on the normalization higher 
than zero. We calculated the errors on the line normaliz¬ 
ations freezing the line energies at the laboratory frame val¬ 
ues, obtaining normalization values significantly lower than 


the values reported in Tableland physically more consistent 
with the expected widths implied by the blurring parameters 
of the fluorescence iron transitions (Table O. 

Finally, we added a zero-width line at ~ 7.11 keV, find¬ 
ing best-fitting position and normalization consistent with 
the values reported in Table O and obtaining a final best-fit 
reduced \ 2 of 1.025 (519 dof) and a satisfactory pattern of 
residuals as shown in Fig. [10] In Fig. llll we show the broad¬ 
band data/model ratio for a best-fit continuum without 
emission/absorption features (panel A); for the best-fitting 
model of Table [3] when the reflection fraction is set to zero 
(panel B); when the depth of the fundamental cyclotron line 
is set to zero (panel C); when the depths of the two harmon¬ 
ics are set to zero (panel D). 

During the fitting procedure we noted that the inner 
disk radius, Rm, of the reflection component becomes rather 
unconstrained by the fit because of a strong correlation 
with the other smearing spectral parameters (i.e. inclina- 
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Table 2. Emission lines in the spectrum of GRO J1744-28. Errors given at Ay 2 = 2.706 (90% confidence for the single interesting 
parameter). 


Ion (Transition 0 ) 

E lab 

E 0 6s 

Width 

Flux'* 

EQW C 


keV 

keV 

eV 

(10- 4 ) 

eV 

Ka Sxvi 

2.62 

2 fi | 5 + 0 - 04 
z.uu_Q 03 

=Caxix 

160l 4 ° 

60l 2 ° 

Ka Ar xvm 

3.32 

0 oq+0.03 

o.zO—0.02 

=Caxix 

150tf 0 

70±5 

Ka Caxix 

3.90 

3.91T0.02 

240±10 

82±{° 

44±6 

Ka Fexxv 

6.70 

6.69T0.04 

210t 2 ° 

67 ±f 5 

55±15 

Ka Fexxvi 

6.97 

6.96T0.07 

=Fe xxv 

40^22 

30±10 

K/3 Fe xxv 

7.88 

7.88 d 

=Fe xxv 

7-9 

7.4T2.5 

K/3 Fexxvi 

8.25 

8.25 d 

=Fe xxv 

5-itI;3 

c 1+2.3 

systematic 


2.264T0.004 

0 d 

60±6 

18±3 

systematic 


7 11 +0-2 

0 d 

12±i 

9.4T2.4 

a Rest-frame energies from Verner et al. 

(19961. 




b Total area of the Gaussian (absolute value), in units of photons/cm 2 /s. 
c Line equivalent width. 
d Frozen parameter. 



Figure 10. Data, unfolded model with additive components (up¬ 
per panel) and residuals in units of a (lower panel) using the 
disk+relxill modell (Table[3]l and a set of ionized emission lines 
with common blurring parameters. 


tion angle and e). Fixing the inclination angle to 18°, 22°, 
and 25°, we obtained a lower limit on Ri n of 90 R 9 , 80 R 9 , 
and 70 R 9 , respectively. For higher values of the inclination 
angle R ln is not constrained. If the value of the e parameter 
is frozen to 3 (a value that would correspond to an isotropic 
point-like irradiating source), we obtained a lower limit for 
Rin of ~ 100 R 9 ; if both the inclination is fixed to 22° (up¬ 
per l imit based on the calculations from iRappaport fc Josa 
1 19971 . Table 3) and e = 3, then the lower limit for Ri„ is 120 
R-g- 

in Table [3l first column, we report the best-fitting res¬ 
ults with Ri„ fixed to a reference value of 100 R 9 . Finally, we 
note that if no systematic error is assumed for the EPIC/pn 
data, then a strict lower limit on Ri„ is derived at ~ 60 
R 9 (Table [3] second column). We obtained similar results 
adopting the EPIC/pn spectrum without the three central 
columns (Table [3] third column). 


Differently from the nthcomp Comptonization model, 
the relxill continuum has no low-energy roll-over, imply¬ 
ing a higher flux of this component in the softest X-ray band. 
This affected the determination of the disk emission, and 
we obtained a shift of the ezdiskbb temperature to lower 
values and a higher value for the disk inner radius (R ez ~ 
90 R 9 ). We found that the reflection component requires a 
high metallicity, where iron over-abundance is ~10 the solar 
value. Also the emission lines from the highly ionized trans¬ 
itions of Si/Ar/Ca were not fully well fitted by the broad¬ 
band reflection model and the residuals strongly required 
more flux at the line energies that we interpreted as an in¬ 
dication of possible higher metal abundance. The comp an- 
ion star of GRO J1744-28 is an evolved giant dGosling et al.1 
l2007l l whose external layers are probably rich in the ejecta 
of the supernova explosion that generated GRO ,11744-28 
and a high level of metallicity in the accretion flow could be 
expected. 

The best-fitting parameters and associated errors are 
shown in first column of Table [3] while in the second and 
third columns of the same table, we comparatively show the 
two best-fitting models using the EPIC/pn data without 
any systematic error (second column), and the EPIC/pn 
spectrum extracted from the wings of the PSF (excising 
the three hottest RAWX columns) and no systematic error 
(third column). The range of best-fitting values and associ¬ 
ated errors is useful to comparatively visualize the impact 
of the assumptions based on the choice to assign a system¬ 
atic error to the EPIC/pn data and the impact of possible 
residual pile-up. 


3.5 The pulsed fraction in the iron region: 
comparison with best-fitting models 

We have shown in Sect. nm that there is a clear drop in 
the pulsed fraction in the iron region of the spectrum (see 
Fig. El) and that the shape of the drop is consistent with a 
b road Gaussian prof ile. Such drop has been also noted in 
iNishiuchi et al.l (Il999l ). who proposed its origin in the not- 
pulsed flux of a relativistically smeared iron line. We have 
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Figure 11. 2-60 keV data/model ratios: for a best-fitting continuum without reflection and gabs (panel A), for the best-fitting continuum 
of Table [3] with Refl. fr. set to 0 (panel B); for the best-fitting continuum of Table l3l with Di set to 0; for the best-fitting continuum of 
Table [3] with D 2 and D 3 set to 0. 


shown that the large iron residuals can be interpreted as re¬ 
flection from a truncated disk and we establish here if the 
shape and strength of the residuals of Fig. [5] are consistent 
with this origin. To this aim, we divided the energy spectral 
range in bins of 100 eV width from 5.5 keV to 7.5 keV. We 
calculate for each channel the observed flux of the reflection 
component, relxill, and of the reflection component ac¬ 
cording to the best-fitting model of Table. [3j We calculated 
the not-pulsed excess flux (F ex ) for the amplitude versus 
energy plot using a best-fitting third-order polynomial func¬ 
tion fit to the pulsed fraction outside the 5.5-7.5 keV range 
as explained in Sect. 12.21 according to the following expres¬ 
sion: 

I'v.x - - l) x Fiat (4) 

where f exp is the value of the expected pulse fraction 
in that energy bin from the best-fitting function, f 0 b s is the 
observed pulsed fraction, and Ftot is the continuum flux per 
bin. In Fig. 1121 we show the results of this procedure. We 


found that the not-pulsed flux excess is in reasonable agree¬ 
ment with the flux of the reflection component. 


4 PULSE-PHASE RESOLVED 
SPECTROSCOPY 

Based on the best time-averaged pulse period during the 
persistent emission, we extracted ten pulse-resolved spec¬ 
tra, equally phase-spaced, to study the spectral evolution 
as a function of the spin phase. We used only data from 
the EPIC/pn instrument in the 2.0-9.5 keV range. To get 
a preliminary idea of the spectral changes, we first plotted 
the count difference between the time-averaged spectrum 
and the phase-selected spectrum as a function of the en¬ 
ergy channel. The result of this operation is shown in the 
heat map of Fig l 131 where no evident substructure is vis¬ 
ible, implying that most of the spectral changes is due to a 
slowly changing slope of the overall continuum. To fit the 
spin-resolved spectra we adopted the best-fitting model of 
Table [3] and we found a global acceptable fit (reduced \ 2 
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Table 3. Best-fitting parameters and errors for models using a self-consistent reflection component to fit the highly ionized emission 
lines. In the first column we assigned to the EPIC/pn channel a systematic error of 0.5%, in the second column the EPIC/pn spectrum 
has no systematic error, in the third column we used the EPIC/pn spectrum obtained by excising the PSF after removing the three 
central columns. Errors given at A \ 2 = 2.706 (90% confidence for the single interesting parameter). 



0.5% sys. 

No sys.error 

3 cols 

N h (10 22 cm -2 ) 

6.4T0.3 

6.15i0.15 

6 2”*” 0 ' 2 
°- z -0.3 

r 

u.oo—o 07 

0.295i0.015 

0.49i0.07 

Norm. a 

0.45i0.06 

0.42i0.01 

u.oo—o.06 

E cut (keV) 

8.56i0.15 

8.48i0.10 

8 94 - *" 011 

o-y^_0.l6 

kT diafc (KeV) 

0.29T0.04 

0.30i0.02 

0.28i0.02 

Rez (Rg) 

QQ+100 

yu -50 

65l®° 

too til 

Ei (keV) 

4.68i0.05 

4.75(t°;“ 

4.59i0.06 

Wi (keV) 

0.68i0.08 

1.02i0.05 

0.84i0.07 

Di (x 1CT 2 ) 

8.7T0.2 

22i3 

16i4 

E 2 (keV) 

10.4i0.10 

ll.4i0.08 

9.43i0.09 

W 2 (keV) 

0.66i0.10 

1.0i0.4 

u - y4 -0.08 

D 2 (x 10" 2 ) 

14i9 

24il5 

25i7 

E 3 (keV) 

15.8±J;? 

16.0i0.07 

15.4i0.7 

W 3 (keV) 

2-e±S:S 

2.3i0.7 

2.6i0.8 

D 3 (x 10- 2 ) 

90i30 

93i22 

95i25 

e (emissivity) 

i c+2.0 

1 1 + 1,2 

o o + l-2 
z '°— 0.8 

Inclination (deg) 

qo+20 

Z °-10 

28i4 

23±4° 

R in (R<?) 

lOO 6 

>60 

100 6 

Relxill Log(£) 

3.09i0.02 

3.09i0.02 

3.1i0.03 

Relxill Ape 

>8 C 

>9 C 

>7 C 

Relxill Refl. fr. (1CT 2 ) 

5.3T0.4 

r O+0.7 

7 fi”*” 1 ' 2 

Sxvi norm/ 

7i5 

5.6il.9 

4.9i3.9 

Ar XVIII norm/ 

8^4 

7.9il.8 

7.3i2.5 

Caxix norm/ 

8l 3 

7.3il.3 

6.2i2.2 

Caxx norm/ 

6i3 

5.4il.3 

7.7i2.0 

Ci PN/RGSol 

1.15i0.02 

1.16i0.01 

1.07i0.02 

C 2 PN/RGSo2 

1.02i0.03 

1.10i0.02 

1.02i0.02 

C 3 PN/JEMX1 

1.19i0.02 

1.18i0.02 

1.36i0.02 

C 4 PN/JEMX2 

1.32i0.02 

1.32i0.02 

1.52i0.02 

C 5 PN/ISGRI 

1.20i0.03 

1.16i0.03 

1.41i0.05 

Xred ( dof ) 

1.025 (519) 

1.408 (519) 

1.284 (519) 


“ In units of photons keV -1 cm' 2 s —1 at 1 keV. 
b Frozen parameter during fitting search. 
c Hard limit for this parameter set to 10. 

d Line normalization in units of 10 -4 photons cm -2 s —1 in the line. 


1287/1241) by setting the following parameters free to vary: 
the photon index (F), the normalization of the relxill com¬ 
ponent, and the reflection fraction. The remaining paramet¬ 
ers of the fit were kept fixed to the best-fitting values of the 
pulse-averaged spectrum (Table QT[l. 

In Fig. 1111 we show the variation of the photon-index, 
reflection fraction and absorbed flux as a function of the spin 
phase, all following the smooth sinusoidal trend of the pulsed 
flux. We clearly observe that the lowest value of the photon- 
index is found in correspondence with the pulse maximum, 
while the highest reflection fraction is at the pulse minimum. 
Because the reflection fraction is the ratio between the reflec¬ 
ted and the direct, primary flux, this result suggests that the 
reflected flux does not sensibly depend on the pulse phase 


and we verified this by calculating the reflected flux for each 
phase-selected spectrum, finding it constant within a ~ 5% 
relative uncertainty. We did not detect statistically signific¬ 
ant variations of the fundamental cyclotron line parameters 
along the pulse phase. 


5 DISCUSSION 

In this paper we presented a complete analysis of the broad¬ 
band (1.3-70 keV) spectral characteristics of the persistent 
emission of the bursting pulsar GRO J1744-28 based on a 
long continuous XMM-Newton observation and on the hard 
X-ray coverage provided by two of the instruments on board 
of INTEGRAL (JEMX and ISGRI). 
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Figure 12. Red line: the reflection component flux according to 
the best-fit model of Table [3] Black data: excess flux from the 
amplitude fraction according to Eq. [4] 


Figure 14. Best-fitting parameters 1 a errors for the ten phase- 
resolved spectra. From top to bottom: relxill 0.5-10 keV un¬ 
absorbed flux in units of 10 —8 erg cm -2 s —1 , photon-index, and 
reflection fraction. 


Contour Map of the Spin-resolved Spectra 
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the strong effective area derivative of the Au-M edge and a 
narrow line at 7.11 keV, where there is a physical spectral 
discontinuity due to the sharp neutral iron edge. The RDPHA 
method assigns a ~ 20 eV uncertainty in the photon energy 
scale reconstruction, and the appearance of such narrow fea¬ 
ture could be caused by a slight mismatch in the energy scale 
reconstruction of photons in this band. 

However, despite these caveats, the overall picture and 
interpretation of the main physical processes producing the 
X-ray spectrum remains quite solid, because it was found 
not to sensibly depend on these instrumental uncertainties 
(Table |3). 


5.1 Continuum emission and cyclotron lines 


Phase 

Figure 13. Heat map for the EPIC/pn spin-resolved spectra. We 
plot on the y-axis the count/s/keV difference with respect to the 
time-averaged spectrum. 

The EPIC/pn spectrum is, however, affected by some 
level of uncertainty, because of the high registered count 
rate causing instrumental spectral distortions and shifts in 
the energy scale (i.e. charge tr ansfer inefficiency and X-ray 
loading, see lPintore et al1l2014Tl and by some level of photon 
pile-up. To process the EPIC/pn data we adopted the RDPHA 
pipeline as suggested by the EPIC/pn hardware team, and 
we obtained our best-fitting results considering both an ex¬ 
cised PSF and the total PSF (full spectrum, see Table [3]) ■ 
The main spectral differences in the two spectra affect the 
determination of the power-law photon index, the cut-off en¬ 
ergy, and the depth of the first harmonic, which resulted each 
other not compatible. Moreover, the EPIC/pn spectrum is 
according to us affected by two other (possible) instrumental 
artefacts: a narrow emission line at ~ 2.2 keV, observed in 
the spect ra of bright sources either in emission o r in ab¬ 
sorption dHiemstra et al.l 1201 ll : IPintore et al .1120141 ). due to 


The X-ray continuum spectrum of an accreting pulsar is gov¬ 
erned by different physical process that are all sensible to the 
local condition of the accreting plasma and to the geometry 
and strength of the NS magnetic field. Soft photons pro¬ 
duced by bremsstrahlung, synchrotron and thermal, black- 
body, emission are Compton up-scattered to higher energies 
either by the energetic, free-falling, electrons (bulk-motion 
Comptonization) or through thermal Comptonization. X- 
ray pulsars accret ing above the critical luminosity ( L cr u 
«3x 10 36 erg s _1 iBasko fc SunvaevHl97fil : iBecker k. Wolfil 
l2005l l have a radiative shock, formed and sustained by the 
Eddington radiation pressure, and a significant part of the 
scattering process is due to thermal Comptonization. For 
sources below the critical luminosity, bulk-motion (or fist- 
order Fermi energization) dominates. Sources in this fat¬ 
ter category show steeper spectra (T >2), whereas sources 
of the former class appear flatter (P < 2, IBecker fc Wolfil 

<2qq3)). 

Because of the complexity and non-linearity of all these 
processes, the pulsar spectra have usually been fitted with 
phenomenological models, that, however, provided in most 
cases a statistical adequate description of the data, and a 
general common frame for comparing the spectra of dif¬ 
ferent sources and detect, even very broad local absorp- 
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tion/emission features dHeindl et al.ll2004h . We found that 
the broadband spectrum GRO J1744-28, however, cannot be 
described in terms of simple power-law models (whatever the 
prescription for the exponential cut-off), or combinations of 
black-body and cut-off power-laws, that are more typical for 
bright un-magnetized accreting NSs. We observed a common 
trend of residuals, extending in the 4.0-10.0 keV range, giv¬ 
ing rise to a spectral curvature that can be satisfactorily 
accounted by a set of broad features in absorption at the 
edges of the iron Ka region, together with a set of moder¬ 
ately broadened features in emission, that can be satisfact¬ 
orily described with a self-consistent reflection component, 
smeared by dynamical and relativistic effects produced at 
the inner edge of a truncated disk. The presence of these 
features is statistically required independently from the ad¬ 
opted continuum. We have modelled it first with the sum of 
a soft thermal disk emission and a thermal Comptonization 
model (nthcomp), that has two main parameters that de¬ 
termine the cut-off, high-energy, roll-over (the electron tem¬ 
perature and the cloud optical depth, related to the T para¬ 
meter). We note here that we obtained very similar results, 
in terms of x 2 and spectral parameters for the local emis¬ 
sion/absorption features, using a power-law with a highecut 
cut-off, a model widely adopted for describing the hard spec¬ 
tra of X-ray pulsars. The addition of an other continuum 
component (as a black-body emission) resulted in thermal 
temperatures that closely followed the shape of the nthcomp 
soft seed-photon temperature, and we concluded that the 
nthcomp encloses all the characteristics of the most com¬ 
monly models used in the literature to fit the spectra of 
bright X-ray pulsars. 

In Sect. l3.4l we showed results from a simpler continuum 
made with a cut-off power-law (relxill model), a softer 
thermal disk emission (ezdiskbb), and a broadband reflec¬ 
tion component. The reflected spectrum is also an important 
continuum contributor, through continuum Thomson scat¬ 
tering and bremsstrahlung emission, and significant changes 
in the best-fitting parameters of the other local features in 
absorption were expected (compare Tableland Table [3j. It 
is clear that the modelling of the underneath continuum is 
more important than the statistical uncertainties from the 
best-fitting modelling, although the general picture of the 
cyclotron features appears quite solid. The uncertainty on 
the cyclotron fundamental energy is of the order of 6%, the 
width of 25%, whereas the depth is the most sensible para¬ 
meter, due to the strong correlation with the general shape 
of the fluorescence iron line. 

If we use the best-fitting parameters from Table [3] for 
the fundamental cyclotron line energy (Ei=4.7 keV), we 
can derive the magnetic field of the NS using Eq.0 leading 
to a magnetic dipole field of (5.27 ± 0.06) X 10 11 G (magnetic 
dipole moment /r = 5.27 x 10 29 G cm 3 ), assuming the CRSF 
is formed at the surface of the NS (gravitational red-shift 
« = 0.3). 

This estimate allows to constrain the magnetospheric 
radius (Rm), where the ram pressure of the accretion flow is 
b alanced by the m agnetic pressure, according to the formula 
dFrank et al.ll2002h : 

Rm = 2.9 x 10 S KmM 1 N / ^Re 2/7 L- 7 2/ G4o 7 cm ( 5 ) 

where K m is a coefficient related to the accretion flow 
geometry (1 for isotropic radial accretion), Mns is the NS 


mass in units of Mq, Re is the NS radius in units of 10 6 cm, 
I/37 is the source luminosity in units of 10 3 ' erg s _ , and /730 
is the magnetic dipole moment in units of 10 3 ° G cm 3 . From 
our B estimate, the bolometric luminosity derived by extra¬ 
polating the observed flux, assuming isotropic emission, and 
classical values for the NS mass and radius (1.4 Mq, 10 6 
cm radius), we derive a value of Rm — 10 8 XK m cm. The 
co-rotation radius for a NS of period P sp in is 


( GMP 2 in \ 1/3 


(6) 


and corresponds for GRO J1744-28 to 10 8 cm. A value 
of Km close to 1 (isotropic accretion) would imply that the 
NS is accreting at the spin equilibrium. However, we derived 
from the best-fitting parameters of the disk emission and of 
the reflection component an inner disk radius lower limit 
closer to a value of 2 x 10 7 cm, that sets the Km lower limit 
value at ~ 0.2, that is in agree ment with the value ~ 0.2 
reported for low B-field pulsars in lBurderi et al.l (l2002l l . and 
not unre alisti cally lower than a more recent estimate of 0.5 
found bvlLong e t al.l (120051 1. We also note that, as discussed 
bv lBozzo et al.l ( 2009l l. the expected value of this parameter 
is affected by large uncertainties due to the different theor¬ 
etical prescriptions for the magnetospheric radius published 
so far. 

We found that the width of the fundamental CSRF is 
15-20% the value of the energy position, and this is in line 
with wh at observed in simi lar other accreting pulsars (see 
Fig. 5 in lHeindl et alJfeoQdl l. Adopting the classical formu¬ 
lation for the ther mal broadening in a plas ma with electron 
temperature kT e llMeszaros fe Nagelll985ll : 



FWHM 


/ 81n2 kT e 
V m e c 2 


1/2 

cos# 


( 7 ) 


assuming an angle 6 ~ 20°, we can derive an electron 
temperature of ~ 11 keV, which is quite close to the cut-off 
energy of the broadband spectrum, estimated at ~ 8.6 keV. 
From the spin-resolved analysis, we found that the funda¬ 
mental cyclotron feature did not show any statistical s igni¬ 
ficant correlation with the pulsed emission. iDeeenaar et al.l 
d2014l l noted in the C/iandra/HETGS observation of GRO 
J1744-28 an unidentihed absorption line at ~5.11 keV that 
was tentatively associated with a cyclotron absorption fea¬ 
ture, though the line width could not be constrained and 
it was kept fixed to 10 eV in the fit. Because the energy 
and the width of the fundamental cyclotron line as derived 
by our analysis of the EPIC/pn data (Table O points to¬ 
wards a significantly broader feature, we retain that these 
two detections are not related. 

We found a suggestive, but statistically weak, hint of 
two cyclotron harmonics at ~ 11 keV and ~ 15 keV. We note 
that lYounes et al.l (l2015h reported also a broad dip in the 
NuSTAR spectrum of this source at ~ 10 keV, whose shape 
resembles the one shown in the INTEGRAL/ JEMX data. 
The widths and depths of these two features could not be 
independently determined, because the features are broad 
and they spectrally overlap. Moreover, we noted a certain 
dependence based on the assumptions for the continuum 
and for the systematics that affect the inter-calibration of 
the different instruments that were used to constrain these 
shapes. Notwithstanding these limitations, the energies of 
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the first and second harmonics were always found close 
to the expected harmonic ratios 1:2 and 1:3 (for a funda¬ 
mental line at ~4.7 keV), and GRO J1744-28 would be, 
after 4U 0115+64 dMuller et al.l 2 0131. CRSF fundamental 
at 14 keV) and V 0332+53 I Pottschmidt et ahll2005l . CRSF 


fundamental at 26 keV), the third X-ray pulsar to show 
a second harmonic in its spectrum. The slight discrep¬ 
ancy in the harmonic ratios (2.2 and 3.3 for the first and 
second harmonic, respectively) has been observed also in 
other X-ray pulsars, and it has been show n to be depen d¬ 
ant on the accretion rate as discussed in IN aka lima et al.l 
(120101) . We found that the two harmonics have larger depths 
with respect to the fund amental, as reported in other 
sources (e.g. for X0115+63. ISantanaelo et al.l 119996 . an ob¬ 
servational Jact_e3gDlained__through the two-photon process 
dAlexander fe Meszarosl 1 199 ill . Because these features he 
outside the EPIC/pn band, we could not track parameter 
variations with respect to the spin-phase. 


5.2 The highly ionized reflection component: 

GRO J1744-28 between the LMXB and the 

HMXB class 

The set of highly ionized feature is satisfactorily compat¬ 
ible with a reflection spectrum from a cold disk truncate at 
large distance from the NS. The inner disk radius estimate 
is correlated with the values of the other smearing paramet¬ 
ers, but for realistic values of both inclination angle (~ 22 
deg) and index of the emissivity law profile (e = 3), a lower 
limit value on Ri n is obtained at 120 R 9 . The temperat¬ 
ure and the luminosity of the disk continuum component 
was f ound com patible with this estimate. Differently from 
[Nishiuchi et alj (Il999h , we do not find evidence of extreme 
red wings, that would imply low inner disk radii, and our 
estimate is rather close to an independent estimate based 
on the iron l ine shape observed in the Chandra observation 
analysed bv lDegenaar et all 1 2014 1. 

The presence of reflection components in highly mag¬ 
netized accreting pulsars has been usually neglected, 
because of the large distance of the inner disk ra- 
dius (if present ) and expected low reflection fraction. 
iBallantvne et all (l2012f) first proposed a set of self-consistent 
reflection tables and showed that moderately broad emission 
features could be signatures of such component as in the case 
of the accreting pulsar LMC X-4. In that case however blur¬ 
ring parameters derived from best-fitting did not provide 
physically consistent results, as the derived inner radii were 
too low compared with the expected disk truncation radius 
at the magnetosphere. On the contrary, for weakly mag¬ 
netized accreting X-ray pulsars the detection of reflection 
components has been already established in a considerable 
amount of cases. A relativistic iron line profile has been first 
reported for SAX J1808.4-3658, where the blu rring paramet¬ 
ers w ere found in line with the expectations (Papitto et al. 
200S), and later also for IGR J17511-3057 i Papitto et al. 
20101 . IGR J 17480-2446 dMiller et al.l l201lll and in HETE 
J1900.1-2455 (IPapitto et al.l l2013l l . In the above-mentioned 
cases the inner disk radii were constrained in the range 20- 
40 R g and were consistent with the possible magneto-spheric 
radii derived from the study of spin timing, whereas in non¬ 
pulsating accreting neutron stars signific antly lower values 
were derived (6-15 R g , see lCackett et alll2010i l. 


We derived an inclination angle between 18° and 48°, 
that is within the range of values for the system accord- 
ing to the possible evo lutionary histories of the system 
llRappaport fe Josa 1997lf. From th e X-ray mass function of 
GRO J1744-28 (jFinger et al.Hl996h . we thus derive an upper 
limit on the companion mass between 0.23 M@ and 0.27 Mq, 
for a 1.4 and 1.8 Mq NS mass, respectively. This range en¬ 
compasses the most probable v alue (0.24 Mq) indi c ated f or 
the present companion mass bv lRappaport fe Jossl dl997l l. 

The reflected component has a high ionization degree 
(log£ ~ 3.1) as expected from the presence of H- and He- 
like emission lines, and from the distance of the inner disk 
radius and bolometric source luminosity, we can estimate the 
reflector electron density n e =Lboi/{r 2 £) ~3 x 10 20 cm -3 . 


6 CONCLUSIONS 

We presented results from the spectral analysis of the broad¬ 
band spectrum of the bursting pulsar GRO J1744-28. Using 
XMM — Newton/INTEGRAL data we covered the 1.3-70 
keV energy range, when the source was accreting near the 
Eddington limit for a 1.4 Mq NS. Independently from the 
adopted continuum, we found in the spectrum a moderately 
broadened absorption feature in the range 4.6-5.3 keV and 
a set of moderately broadened emission features. We inter¬ 
preted this absorption dip in terms of a cyclotron resonant 
scattering feature and we also found suggestive evidence of 
two other absorption features close to the energies expected 
for the first and second harmonic. From this interpretation, 
we derived an estimation of the GRO J1744-28 pulsar mag¬ 
netic field, 5.27 + 0.06 x 10 11 G. 

We located the region responsible for the moderately 
broadened emission lines outside the pulsar magnetosphere, 
in a reflecting, accretion disk. The study of the pulse amp¬ 
litude fraction and the spin-resolved spectra indicated that 
the iron complex in emission is not pulsed and its flux is 
not phase-dependent. We found that the disk surface is 
highly photo-ionized, with log(£) ~ 3.1, with a metal over¬ 
abundance of a factor ~ 10 and truncated at a distance of 
> 100 R g , where dynamical effects, caused by the high velo¬ 
cities of matter in the disk, are still important. The spectrum 
is also characterized by soft disk emission, whose trunca¬ 
tion radius results compatible with the estimate from the 
reflected component. These two independent, but conver¬ 
ging, constraints allowed us to estimate the Alfven efficiency 
parameter for disk-accretion, close to 0.2. The lower limit on 
the inclination angle of the disk is 18° and compatible with 
independent e stimates on the possible evolutionary history 
of the system (IRappaport fe Josslll997ii . 
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